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METHOD FOR CALCUIATION OF COMPRESS IBIE LAMINAR 
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PRESSURE GRADIENT WITH ZERO 
HEAT TRANSFER 

By Morris Mo r due how and Joseph H. Clarke 
SUMMARY 


The Karman-Pohlhausen method is extended primarily to sixth-degree 
velocity profiles for determining the characteristics of the compressible 
laminar boundary layer over an adiabatic wall in the presence of an axial 
pressure gradient. It is assumed that the Prandtl number is unity and 
that the coefficient of viscosity varies linearly with the temperature. 

A general approximate solution which permits a rapid determination of 
the boundary-layer characteristics for any given free -stream Mach number 
and given velocity distribution at the outer edge of the boundary layer 
is obtained. Numerical examples indicate that this solution will in 
practice lead to results of satisfactory accuracy, including the critical 
Reynolds number for stability. For the special purpose of calculating 
the location of the separation point in an adverse pressure gradient, a 
short and simple method, based on the use of a seventh-degree velocity 
profile, is derived. The numerical example given here indicates that 
this method should in practice lead to sufficiently accurate results. 

For the special case of flow near a forward stagnation point it is shown 
that the Karman-Pohlhausen method with the usual fourth-degree profiles 
leads to results of adequate accuracy, even for the critical Reynolds 
number. 


INTRODUCTION 


In reference 1 it was concluded that from the viewpoint of both 
accuracy and convenience of calculation a suitable method for deter- 
mining, the characteristics of a compressible laminar boundary layer is 
that based on an extension of the Karman-Pohlhausen integral method to 
velocity profiles of higher degree than the fourth, especially sixth 
degree. An ordinary differential equation for general types of flow was 
derived, but only the flow over a flat plate at zero incidence was 
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investigated in detail. The purpose of the present investigation is to 
apply explicitly this method to flows with axial pressure gradients. 

The ordinary differential equation derived in reference 1 is con- 
verted here into a convenient nondimensional form, and a general approxi- 
mate solution of this equation in closed form is then derived. By means 
of this solution the physically significant- boundary-layer characteristics 
of the flow over an insulated wall can be calculated fairly easily and 
quickly for any given free -stream Mach number and for any given velocity 
distribution outside of the boundary layer. For particular accuracy in 
determination of the separation point in an adverse pressure gradient, 
a simple method based on the use of a seventh- degree velocity profile, 
which, in accordance with a suggestion of Timman (reference 2), is made 
to satisfy an additional condition at the separation point, is derived. 

For illustrative purposes, two simple but basically different types 
of flow are treated in detail. In the first example flow with a linearly 
diminishing velocity at the outer edge of the boundary layer is con- 
sidered. Such a flow is of particular interest here because it repre- 
sents the simplest case of an adverse pressure gradient and because the 
usual application of the Karm^n-Pohlhausen method with fourth- degree 
velocity profiles has been known to lead to., highly inaccurate results 
in such cases (cf., e.g., reference 3)* The skin friction, velocity 
profiles, and critical Reynolds number are calculated for several Mach 
numbers in this case by means of the general solution derived here. 

These results are then compared with known series solutions of the 
partial differential equations (references 4 and 5) for this case; the 
agreement is shown to be satisfactory for practical purposes. This 
agreement includes critical Reynolds numbers for l8minar-flow stability 
and location of the separation point. In connection with the latter it 
is significant to note that for zero Mach number the Karman-Pohlhausen 
method with fourth-degree velocity profiles leads to an error of about 
30 percent in predicting the separation point, while the use of a sixth- 
degree profile reduces this error to 15 percent, and the special use of 
a seventh-degree profile for this (and only this) purpose practically 
eliminates the error for all Mach numbers (0 to 3) considered numeri- 
cally here. 

The second example treated is the flow near a forward stagnation 
point. This case is of interest not only because it may represent the 
subsonic flow over a blunt nose but also because it had already been 
found by Schlichting and Ulrich (reference 6) that in such a case no 
physically significant values can be found for the boundary- layer thick- 
ness when sixth-degree profiles are used. It is shown in the present _ 
investigation that although this is, strictly speaking, true, an approxi- 
mate value for the thickness is nevertheless obtainable. It is also 
shown that in this case the usual type of application of the Karman- 
Pohlhausen method' with fourth-degree velocity profiles leads not only 
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to real values of the boundary-layer thickness but also to results of 
satisfactory accuracy, even for the critical Reynolds number. Conse- 
quently, in such a case fourth-degree profiles may be used with confi- 
dence in the accuracy of the results to be obtained. 

As in reference 1, it has been assumed that the viscosity coeffi- 
cient is directly proportional to the absolute temperature, and that, 
following a suggestion of Chapman (references 7 and 8), the constant of 
proportionality may be chosen so that Sutherland's relation is exactly 
satisfied at the wall. This introduces considerable mathematical sim- 
plification and will lead in practice to fairly accurate results at 
least for lower Mach numbers (below 5)* For the same reasons it has 
also been assumed that the Prandtl number is unity, and this may be 
viewed as an approximation for air, where the Prandtl number is more 
nearly 0 . 72 . 

This work, carried out at the Polytechnic Institute of Brooklyn 
Aeronautical Laboratories, was sponsored by and conducted with the 
financial assistance of the National Advisory Committee for Aeronautics. 


SYMBOLS 


a 2 

a 2 

b 

C 

C f 


Fl,F 2 > F 3 


p 

coefficient of t in velocity profile; also given 
by equation (20) 

constant "average" value of a 2 

constant in given velocity distribution u p > 

equations ( 27 a) and (30a) 

factor of proportionality in equation = c(t/ T^) 

skin-friction coefficient (*• du/dy) D (Woo 2 / 2 )) 

specific heats at constant pressure and at constant 
volume, respectively 

integrals defined by equations (12) 

values of F-]_, F 2 , and F^ based on seventh-degree 

profiles; given by equations (B2) 


■^1s?^2s^ f 3s 
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constant defining strength of a shock wave in flow 
outside of boundary layer; given by equation (l9e) 

integrating factor of equation (21) 

constant defined by equations (l8a) and (l8b) 

thermal conductivity of fluid 

characteristic length; see equations (13) 

characteristic length in stability calculations 
(appendix C) 

Mach number 

pressure 

Reynolds number 

free-stream Reynolds number (po, 

Sutherland's constant, defined by equation (8); for 
air, S = 216° R 

absolute temperature 

transformation variable defined by equation (10) 

velocity components in x- and y-directions, respectively 

coordinates parallel and perpendicular to surface, 
respectively 

quantities defined by equations (C2) 
ratio of specific heats ( G p/ C v)j *’ or 7 = 1.4 

physical boundary- layer thickness in xy-plane 
boundary-layer thickness in xt-plane 
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Subscripts 


value of A based on seventh-degree profiles and used 
only for determining separation po int 

coefficient of viscosity 

kinematic viscosity (p/p) 

dimensionless distance along wall (x/l) 

mass density 

dimensionless variable (t/S-^) 


0 values at wall, for example, T 0 

1 local values at outer edge of boundary layer; for 

example, M]_, Tj_, and Uj_ 

a initial value in region of adverse pressure gradient 

(appendix B) 

b values immediately behind a shock wave at leading edge 

c values at critical point for stability (appendix C) 

s value based on seventh-degree profiles 

<» values in undisturbed free stream; for example, 

H», and u^ 

A prime ( ' ) denotes at first differentiation with respect to x. 
In equation (l4) and thereafter in the main text, a prime denotes dif- 
ferentiation with respect to £. In appendix C, it denotes differen- 
tiation with respect to y/l. 


BASIC EQUATIONS 

Two-Dimensional Compressible Flow 


The following equations characterize the steady two-dimensional 
compressible flow of a gas in the laminar boundary layer over a surface 
of large radius of curvature compared with the boundary- layer thickness: 


8u , „ 8u 
pu ^ + pv T- = 
ox oy 


8u , , 8 ( cJu\ 

5y dy^ dy/ 
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St 

puc * s + 


St _ 

pvc p ^ = 


- p i u i u i' u + 


±L^\ 

dy\ Sy/ 



(2) 


S 

Sx 


(pu) + 



= 0 


( 3 ) 


Pl T 


(4) 


Equations (l) and (2) are the momentum and. energy equations, respectively. 
The equation of continuity is given hy expression (3), while equation (4) 
follows from the ideal gas law. In equations (l), (2), and (4) account 
has been taken of the well-known implication of Prandtl's basic boundary- 
layer assumption that the pressure gradient Sp/Sy across the boundary- 
layer thickness is zero. Thus the axial pressure gradient can be con- 
sidered as given by the Bernoulli equation 


Sp , 

5 = -W 


(5) 


where the subscript 1 refers to the potential flow over the surface 
streamline outside of the boundary layer. This potential flow is 
assumed to be known. 

For a Prandtl number pcp/k of unity and zero heat transfer at the 
wall (8T/8y =0 at y = 0 } it can be shown, by eliminating p^u^u^’ 
from equations (l) and (2), that 


u2 

T 


+ CpT . = Constant 


( 6 ) 


Equation (6) can be expressed in the following nondimens ional form: 


1 + 
T X 


7-1 




2 


(6a) 
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Thus the energy equation (2) can for this case be replaced by equa- 
tion (6a ) , which gives the temperature explicitly as a simple function 
of the velocity. This case of unit Prandtl number and zero heat trans- 
fer is the one which will be treated in the present analysis. 

It will be assumed here, as in references 1, 7> and 8, that the 
viscosity coefficient can be expressed as a linear function of the 
temperature in the form 



(?) 


where C is a constant determined so that relation (7) will give the 
same value of p at the wall as the Sutherland viscosity-temperature 
relation. The latter relation is generally assumed to be accurate at 
least below the hypersonic range. Thus 


e = i * Eft.) . 

VW (t 0 /t«,) + (s/tJ 

where- S is a constant which for air has the value S = 2l6° R. From 
equation (6a) it follows that the wall temperature T 0 as a function 
of the free-stream Mach number M^ will be: 


T 

= 1 + 
T„ 



(9) 


Hence the constant C will be a function of the Mach number M„, and 

00 

the free-stream temperature T^. 

By integrating equation (l) with respect to y across the boundary- 
layer thickness y = 0 to y = 5 and introducing the Dorodnitzyn vari- 
able t defined by the transformation 


y = 




dt 


( 10 ) 
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the following equation, essentially the Karman integral momentum equa- 
tion for compressible laminar boundary layers, can be derived (cf. ref- 
erence 1, equation (24)): 



p iV 


■)' - F 2 (l 0ge 


S 2 = F„ 


Ti 

p l U l T ~ 


( 11 ) 

Here 



( 12 ) 


and S t , the value of t when y = 5, is the boundary- layer thickness 
in the xt -plane,., while t = t/6-fc. 

Introducing the dimensionless variables X and | defined by 



L 


where L is a characteristic length and R M = p^u^l/p^ is the free- 
stream Reynolds number, equation (ll) can be written as 


Pi' 


+ hl/ 2 

f 2 \ 

Pi 

y i 

"i l 

F J 


^ £3 

U 1 Pi F-l 


(14) 


where the primes henceforth denote differentiation with respect to g 
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Velocity Profiles 

If a definite form of profile for u/u^ as a function of t is 

assumed, then equation (l4) becomes an ordinary differential equation 
in In the present analysis, sixth-degree profiles will be 

assumed for most cases and will be chosen to satisfy appropriate bound- 
ary conditions at the wall and at the outer edge of the boundary layer. 
This is in accordance with the results of reference 1, where it was 
concluded that, with the possible exceptions of determination of sep- 
aration point in an adverse pressure gradient and of calculation of the 
boundary layer in the vicinity of a stagnation point, the use of sixth- 
degree profiles should for practical purposes lead to sufficiently accu- 
rate determination of laminar-boundary-layer characteristics, including 
critical Reynolds number for stability. Sixth-degree profiles have also 
been applied recently by Weil (reference 9). 


The following boundary conditions, expressed in the xt-plane, will 
be satisfied: 


At t = 0: 


At t = 1: 


JL = o 

U 1 

d 2 ( u / u l) 

&r 2 


PlM' To ^ 
It-i ) p lu J T 

\ x/ 00 \ 00 / 00 


d 3 (u/ui) 

St 3 


= 0 


(15) 


— = 1 
u l 

d(u/ u i) _ ^ 2 (u/u 1 ) _ 8 3 (u/u-l) 
8t . St 3 


It should be noted that the second condition at the wall (t = 0) in equa- 
tions ( 15 ) follows in general from the original partial differential equa- 
tion (l) in conjunction with equation (7). The third condition at the 
wall follows from differentiation of equation (l) with respect to t, in 

conjunction with the condition of zero heat transfer ((Sr/dT) 0 = o) at the 
wall. The second and third conditions of equations (15), moreover, are 
valid only for zero normal velocity at the wall. 
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The sixth-degree polynomial satisfying conditions (15) is found 

to he 

JL _ 2^1 - + a 2 T 2 - (5 + 2a 2 ) T ^ + 2 (3 + * 2 )^ " ( 2 + f a 2^ T ^ 

where : 


_ iV * o Pi W 2 
2 ’ 2 T to pjpj 


(17) 


Substitution of equation (l6) into equations ( 12) leads to the 
following explicit expressions for F]_, F 2 , and Fg in terms of a 2 : 


F-l = 0.1093 + 0.00211a 2 - 0.000622a 2 2 

F 2 = -0.2857 - 0.01905a 2 " y ~~2^ m i 2 (°-3950 + 
0.021l6a 2 - 0.000622a 2 2 ) 

Fg = 2 - 0.4000a 2 


(18) 


By substituting expressions (l8) and (17) into equation (l4), the latter 
becomes an ordinary differential equation of the first order in Ml)* 

The quantities with subscript 1 may be regarded as given functions 
of | from the potential flow outside of the boundary layer. From the 

relation ( u 2 /2 ) + c p T = Constant, which is valid outside of the boundary 
layer even across a shock wave, it follows that 





(19a) 
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Moreover, the local Mach number M-j_ can be expressed by: 



(19b) 


If there are no shock waves, the flow outside of the boundary layer may 
be considered as isentropic, so that the relation 


Pi 


Ti\l/(7-l) 
T~ ) 

. CM / 


(19c) 


is valid. In the presence of a shock wave, the flow may be considered 
as isentropic before and behind, though not across, the wave. In that 
case, equation (19c) may be modified to 


Pi 



,l/(7-l) 


where G is a constant given by 


(I9d) 


G = 


p b V^V 1/(7 " 1) 

Pm/\1oo / 


(I9e) 


while pb and T^ are 'the values of p and T immediately behind the 
wave, and p-j_ and T-j. refer to the region behind the wave at the outer 

edge of the boundary layer. Relation (19c) can be considered as a special 
case of relation (l9d), namely the case of G = 1. The ratios Pb/P^, and 

Tb/Too may be considered as known from the given flow outside of the bound- 
ary layer. 


a 2 


By the use of equations (9), (19a), and (l9d) the expression (17) for 
can be written, more explicitly, as 


-5 S 1 4 + ^ 1 + (^Kf - ^ 


,3-27 

7-1 


2 


( 20 ) 
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GENERAL SOLUTIONS OF EQUATIONS 


For a given potential- flow velocity distribution uqJu^C |) and 
Mach number M , the differential equation (l4) can be solved for X(£) 

without any basic difficulties hy well-known numerical methods, such 
as Adam’s method or the Runge-Kutta method (cf. , e.g., reference 10). 
Such a straightforward procedure may in practice nevertheless be tedious. 
A relatively simple general approximate solution of equation (l4) will 
therefore be derived in this section. This solution will be found to 
be sufficiently accurate for most practical purposes. 


General Approximate Solution 

The solution to be derived here will be based on the approximating 
assumption that in the expressions for Fp and F£ the quantities a£ 
and a 2 2 as given by equation (20) may be replaced by constant "average" 

values a 2 and a 2 2 for the flow. This assumption is justified by the 
fact that in expressions (l8) for F]_ and F 2 the a 2 terms are rela- 
tively small; hence even large errors in the evaluation of a 2 will lead 
to only small errors in the evaluation of F-j_ and F 2 . 

By replacing a 2 and a 2 2 by constant average values a 2 and 

ap 2 in expressions (l8) for F;l and the quantity F^ becomes a 

constant, while F 2 can be written as 

( 18 a) 


(l8b) 


F 2 = -0.2857 - 0.01905a 2 - ^ 7 - - jKM-) 2 

where K is a constant given by 


K = 0.3950 + 0.021l6ff 2 _- 0.000622a 2 ‘ 
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By use of the relation 



1 + 



which follows from equation ( 6 a ) , and hy insertion of expression (l 8 a) 
for F 2 , expression (l 8 ) for F 3 , and expression ( 17 ) for a 2 in the 

right side of equation (l4), the latter equation can he written in the 
form: 


A* + 2X1 


pil + v 




+ ^0.08571 + 0.01905a 2 ) + 



F 1 u x P] _ T. 


( 21 ) 


Equation (21) is now a linear differential equation of the first 
order in X. By use of relations (19a), (19b), and (l9d) it will he 
found that an integrating factor I of equation ( 2 l) is: 


I 



~ (0 . 08571+0 . 01905a 2 )/ Ti n 


2 

7-1 





The solution of equation (21) satisfying the boundary condition X = 0 
or some finite value (in the case of a stagnation point at £ = 0 ) at 
1=0 is then: 


,^±( 0 . 085 n + 0 . 0190 ^ ( — (0.1950+0. 0211 6a 2 -0.000622s^ 

x _ _4_ Jo iv) \ T coJ 

FlG ^ ^ + JL(° . 08571+0. 0190552^^ - ^-(0.1950+0. 021l6a 2 -0. 000622^5) 

where T^T^ is obtained from equation ( 19 a). 


) 


dl 


( 22 ) 
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For a given velocity distribution u^Ju the effect of Mach 

number is given by the T^Jt^ terms in equation (22) and, in the 

case of a shock -wave at the leading edge, also by the constant G 
(cf. equation (l9e)). If the flow over a given object is considered, 
however, then the Mach number effect will be contained also in the 
ujjugo terms, since the velocity distribution u^j u^ will then be a 

function of Mach number. 

The integral in equation (22) can in practice be evaluated without 
difficulty by numerical means such as Simpson's rule (cf., for example, 
reference 10). An average value a^ for. can, in any given case, 

be chosen by considering equation (20). 


Boundary-Layer Characteristics 


Once X,(i) has been determined, the boundary-layer characteristics 

follow from the equations developed here. Thus the local skin- friction 

coefficient Of will be;l 


(p Bu/Sy ) 0 
C f = - = 4\|C 


P^co 2 / 2 


% Tl j 1 - j) -1/2 


U 00 Too 


(23) 


The physical boundary-layer thickness in the xy-plane can be determined 
from transformation (10), and an explicit expression for it in terms of 
\ is given in appendix A. 

The velocity profiles in the xt-plane are given by equation (l6). 
These can be converted to the physical xy-plane by means of transforma- 
tion (10). An explicit expression for y in terms of t is given in 
appendix A. The temperature profiles follow from the velocity profiles 
by means of equation (6a). It should be noted, however, that for most 
practical purposes it is not actually necessary to determine the profiles 
in the xy-plane* Skin friction (cf. equation (23)) and separation point, 

■^It is demonstrated subsequently that the flow near a stagnation 
point is better represented by a fourth-degree velocity profile than by 
a sixth-degree profile. In such a region the expression for the skin- 

friction coefficient is altered in that the term ^1 - -~j is replaced 

a 2 \ 
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for example, can be determined directly from results in the xt-plane, 
while stability calculations (cf. appendix C) can be performed without 
use of the profiles in the xy-plane. 

The equations developed thus far should be adequate for the prac- 
tical determination of boundary- layer characteristics in most cases. 
There are, however, two important exceptions. These are: (a) The loca- 

tion of the separation point in an adverse pressure gradient and (b) the 
flow in the vicinity of a forward stagnation point. Item (a) is treated 
in the immediately succeeding paragraph, while item (b) is investigated 
in a later separate section. 


Determination of Separation Point 

Although the equations thus far developed will usually be found to 
lead to sufficiently accurate results for most boundary- layer character- 
istics even in flows with adverse pressure gradients, it will be found 
that the location of separation points will still be predicted with 
appreciable error. This error can be greatly diminished by a simple 
modification of the preceding equations. 

For the purpose of determining the separation point, velocity pro- 
files can be chosen which, in addition to satisfying boundary condi- 
tions. ( 15), satisfy a further boundary condition at the separation point . 
The additional boundary condition is obtained by differentiating equa- 
tion (l) twice with respect to t, and by then taking values at the wall 
at the point where du/St = 0. In this manner, with the use of equa- 
tion (6a), the condition 



is obtained. 

Condition (24) is strictly valid only at the separation point (where 
Su/dy = 0). Nevertheless, by using this condition for the entire flow, 

It seems plausible that the location of the separation point in any given 
case could be predicted with improved accuracy, since the value of the 
dimensionless boundary- layer thickness X may be expected to be thereby 
more accurately determined at the separation point. This has already 
been shown by an example for incompressible flow In reference 2 and will 
be further shown by examples (table I) for compressible flow in the present 
analysis. 
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A seventh-degree polynomial in t satisfying conditions (15) and 
(24) can be chosen. By proceeding in the same manner as (previously) 
with the sixth-degree profiles, the following solution \ B for X(|), 

analogous to equation (22), is obtained (a detailed derivation is given 
in appendix B): 


32.83 J 

[ s | 

fa.'* 

6.128| 

i T 

V 00 / 

>0.9366 
| d| 

° ( 

^y.128/^ 

2.437 

1 


(25) 


The value of y = 1.4 has been used in equation (25). In case the 
region of adverse pressure gradient follows a region of favorable pres- 
sure gradient, and hence starts at some point behind the leading eige, 
then equation (25) should "be modified to equation (b 4) in appendix B. 
Separation will occur at the point where (see appendix B): 



The location of the separation point is thus determined at the 
value of | for which equations (25) (or (b 4)) and (26) give the 
same value of \g. For any given adverse pressure gradient (negative 

uj'(|)) and free-stream Mach number M , the function \ s (|) can be 

plotted against | in accordance with both equation (25) (or equa- 
tion (B4) ) and equation (26) in the (anticipated) vicinity of the 
separation point. The point of intersection of these two curves is 
then the separation point. ^ It should be emphasized that equation (25) 


^his point will, according to equations (25) and (26), be inde- 
pendent of G; hence one may put G = 1 for this purpose. This does 
not imply that the separation point will actually be unaffected by a 
shock wave at the leading edge, since the values of u-^u |) and 

tii'/ uj |) may still contain the effect of such a wave. 
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is to be used only for determination of separation point. For other 
purposes, such as determination of skin-friction distribution or of 
critical Reynolds number, equation (22) and its related equations are 
to be used, even for a region in an adverse pressure gradient, unless 
properties close to or at the separation point are sought. 

An application of the general solutions developed here, as well as 
a check on the accuracy of results obtained from them, will now be given 
by means of examples based on two different types of flow outside of the 
boundary layer. 


FLOW WITH LINEARLY DIMINISHING- VELOCITY 


The case 



u 


00 


b| 


(27a) 


where b is a positive constant, is now considered. This represents 
the simplest type of adverse pressure gradient. By introducing the 
linear change of variables 


t 



> 


X]_ = bA 


( 28 ) 


equation (l4) remains unchanged, except that £ is replaced by £■]_ 
and X, by Xp. Moreover, ujju^ now becomes 

— - 1 - !i (27b) 


Thus the constant b no longer appears in the equations. Consequently, 
in this case it is actually necessary to solve equation (l4) only for 
b = 1, since for any other value of b it is necessary only to replace £ 
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by I ]_ and X by X-]_ in the solution obtained for b = 1. The 
present example , therefore, will treat the case 


— = 1 - ! (27c) 

A constant average value a 2 for a 2 can he chosen by noting that 
in this case ag will vary from 0 (at | = 0) to 5 (at the separation 
point, where, according to equation (l6), 1 - (a 2 / 5) = o). Consequently 

the value sT 2 = 2.5 may be chosen. Moreover, the value a 2 ^ = 12.5 will 

be used here. With these values, the approximate solution (22) for the 
case represented by equation (27c) without a shock wave (G = l) at the 
leading edge becomes: 


„,j 

Pi 

^ 5.50 

I ♦ r ~ 2 1 H. 2 U - n 2 ) 

5 / S, 

■^ 6.50 

l + 

( 7 2 X 2 4 ' 

2.75 


(29a) 


where tj = 1 - £. For zero Mach number, that is, = 0, equation (29a) 
reduces to: 


x = 5.76 (1 - 0" 6 ‘ 5 ° - i] 


(29b) 


The skin friction follows from equations (23) and (20), while the 
velocity profiles can be obtained, as previously explained, by means 
of equations (l6) and (A2) (appendix A). For = 1 and 3 , the inte- 
gration in equation (29a) was performed by expanding the integrand in 
a series according to the binomial theorem. The results for X(|), 
skin friction at the wall, and velocity profiles are shown in figures 1, 
2, and 3, respectively, with M m = 0, 1, and 3. 
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The critical Reynolds number for the stability of the laminar 
boundary layer in the presence of small disturbances was calculated 
here on the basis of criteria developed by Lin and Lees (references 11 
and 12). The method of calculation is outlined in appendix C, and the 
results for the critical Reynolds number P 00 u oo^/l J, oo a "^ station 

| = 0.0496 are given in table II for M^ = 0, 1, and 3. By comparing 
the results for zero Mach number with those for the critical Reynolds 
number (1,500,000) at the same station for the flow over a flat plate 
based on the solutions of reference 8 (see also reference l), it is 
seen that the adverse pressure gradient is here highly destabilizing. 
Moreover, it can also be seen from table II that higher Mach numbers 
greatly destabilize the flow over an adiabatic wall. 

As a check on the accuracy of approximate solution (22) to differ- 
ential equation (l4), equation (l4) was solved exactly for this case 
for 1^ = 0 and numerically by Adam's parabolic method (reference 13) 
with increments of = 0.01 for ^ 0. The solutions thus obtained 

for X(£) and for the profiles are shown in figures 1 and 3 , respec- 
tively. The agreement between these solutions and those based on equa- 
tion (22) is seen from these figures to be satisfactory for practical 
purposes. As a further check, the critical Reynolds number at | = 0.0496 
was calculated for M^, = 0 by using the exact solution of equation (l4). 
Table II indicates satisfactory agreement between this result 
(( R co)c = 380,00o) and that ((RjJq = 330, 000 ) based on equation (22). 

It is thus Indicated that equation (22) is a satisfactory approxi- 
mation to an exact solution of differential equation (l4). It is also 
desirable to check whether an exact solution of equation (l4), or 
approximation (22), is a satisfactory approximation to an exact solu- 
tion of the original partial differential equation (l). This would con- 
stitute a check on the practical reliability of the basic method used 
here, namely the extension of the Karman-Pohlhausen method to sixth- 
degree profiles. Such a check can be made by comparing the results 
obtained by means of equations (l4) and (22) with the results of Howarth 
(reference 4) for zero Mach number, based on an accurate solution of the 
original partial differential equations (l) and (3). Such a comparison 
Is shown in figures 2 and 3 for skin friction and velocity profiles, 
respectively, and In table II for critical Reynolds number. In all 
cases the comparison indicates that the solutions obtained by means of 
the equations developed here give sufficiently reliable results for 
practical purposes. 

The separation point is determined by equating the right sides of 
equations (25) and (26). In the present case the point at which these 
right sides are equal was determined by plotting X s (|) according to 
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both equations (25) and (26) on the same sheet in the vicinity of 
| = 0.l4 for = 0 and 1 and £ = 0.10 for = 3. These are the 

vicinities in which, according to figure 1 (based on sixth-degree, not 
seventh-degree, profiles), separation might be expected to occur. The 
separation points determined by this means for = 0, 1, and 3 are 
shown in table I, together with the results of Stewartson (reference 5). 
The agreement between these two sets of results is seen to be excep- 
tionally good. 

For purposes of comparison, the separation point was also deter- 
mined by thb use of equation (22) and by means of the classical Karmln- 
Pohlhausen method based on fourth-degree profiles, the latter for zero 
Mach number only. The results are shown in table I. Although the use 
of the sixth-degree profiles (equation (22)) gives a more accurate 
location of the separation point than the fourth-degree profiles , the 
special use of the seventh degree (equations (25) and (26)) for this 
purpose gives here virtually the exact values. 

It may be asked whether the seventh-, instead of the sixth-degree 
velocity profile could be used for determining the other boundary-layer 
characteristics, in addition to the separation point. Consequently, the 
skin friction and velocity profiles for the present case were determined 
on the basis of the seventh- degree velocity profile, in conjunction with 
equation (l4) (cf. appendix B). It was found that the results, even at 
a station fairly close to the separation point, did not agree so closely 
with Howarth's solution (reference 4) ae those based on the sixth-degree 
profile. It is therefore concluded that the seventh-degree profile 
should be used only as a means of determining the separation point in 
regions of adverse pressure gradients and that otherwise (except for 
flow near a stagnation point) the equations developed on the basis of 
the sixth-degree velocity profile should be used. 


FLOW NEAR A FORWARD STAGNATION POINT 


The case 


— = b£ (30a) 

Uoo 


represents physically the flow in the vicinity of a forward stagnation 
point, such as over the leading edge of a blunt object in subsonic flow. 
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As in the preceding case (equation (27a)) the constant b here can be 
eliminated by means of transformations (28). Hence only the case 


u„ 


(30b) 


need be treated. 

For zero Mach number (m,,, = o) an exact solution of ordinary differ- 
ential equation (l4) in the case represented by equation (30b) is of the 
form X = Constant. This can also be seen from solution (22) for this 

case, with a 2 = a 2 and a rfi = a where according to equation (20) 
(with Me = 0, G = 1, and Uj'/u^ = 1 ) &2 = ~V 2 * From equation (22), 
in fact, the following cubic in X is obtained: 


f(x) = X 3 + 37. 4-3 X 2 - 979. 4x + 6435 = o 


(3D 


By plotting f(X) against X, it is found that there are no physi- 
cally significant roots of equation (31). Such a result was already 
obtained by Schlichting and Ulrich (reference 6). It will be found, 
however, that the curve f(X) against X in the physically significant 
region comes relatively close to the X-axis at a point when the curve 
has a local maximum value. If the value of X at this point is taken 
for purposes of an approximation, then 


X = 9.481 


(32) 


Since the approximation (32) does not appear to be satisfying in 
principle, the use of fourth-degree velocity profiles will be inves- 
tigated for this case. The fourth-degree polynomial satisfying all the 
boundary conditions in equations (15) with the exception of the two con- 


ditions involving the third derivative of velocity 
at t = 0 and at t = 1*) is 


(d 3 (u/ui) / 



= 0 
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With the profile (33) , equation (l4) remains valid, hut the F's aB 
defined hy equations (12) are now given explicitly by the following 
expressions: 


F 1 = 0.1175 + 0.00212a 2 - O.O0O44la 2 2 

f 2 = -0.3000 - o.oi667a 2 - 7 ~ 1 M x 2 (o.4i75 + 
0.018 78a 2 - 0.00044la 2 2 ) 


(34) 


F 3 = 2 - 0.3333s 2 


Expressions (34) replace, in this case, expressions (l8) based on sixth- 
degree profiles. For the case (30b) with = 0, equation (l4) leads 

to a cubic in X whose physically significant root is now found to be 
(cf. also reference 3) 


X = 7.052 (35) 


To compare the accuracy of solutidn (32), based on sixth-degree 
velocity profiles, with that of solution (35); the skin friction, veloc- 
ity profiles, and critical Reynolds number for laminar-flow stability 
have been calculated on the basis of both of the solutions. The results 
thus obtained, together with a comparison with results of an exact solu- 
tion (reference 1.4) for this case, are shown in tables III and IV and in 
figure 4. From these comparisons it is seen that the solution based on 
fourth-degree profiles leads in this case to more accurate results than 
that based, on sixth-degree profiles for the boundary- layer character- 
istics considered, including the critical Reynolds number. It may there- 
fore be concluded that for a stagnation-flow case represented by equa- 
tions (30a) and (30b), a satisfactory solution for practical purposes 
can be obtained by means of equations (33); (34), and (35) based on 
fourth-degree velocity profiles. 
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For M ro £ 0, an approximate solution of equation (l4) quite similar 

to equation (22) can be obtained for the present case. Using expres- 
sions (34) and putting u^u^ = £ in accordance with equation ( 30 b), 

the following approximate solution is thus obtained: 


F-.G 


X = 


P 5 3+ ^-(o.l333+o.oi667ar 2 )/ Tl \^- jr-(o. 2508+0.0x87852-0.000441^) 


£ 


UO 


<15 


4+ ^-(0. 1333+0. 0l667a 2 ) f^Jjz - ^ (0 . 2508+0 . 0l878a 2 -0 . 000441^ ) 


(36) 


where F^ is now given by equation (34). At the leading edge (| = 0 ) 
the value Xq of X will be, according to equation ( 36 ), 




7 - 1 
2 


7-2 


(0.1333 + 0.0l667a 2 + 2 F-l) 


(36a) 


As a numerical example, the case M = 1 was calculated by means 

of equations ( 36 ) and (36a), assuming a 2 = -4.48 and a 2 1 2 = 20.1. The 

results for 4(1), skin friction, and velocity profiles are shown in fig- 
ures 5, 6, and 7, together with the results for = 0. 


CONCLUSIONS 


From the foregoing analysis, based on an extension of the Karm^n- 
Pohlhausen method primarily to sixth-degree velocity profiles for deter- 
mining laminar-boundary- layer characteristics in the compressible flow 
over an adiabatic wall in the presence of an axial pressure gradient, 
the following main conclusions can be drawn. 

1. For a given free-stream Mach number and a given velocity distri- 
bution outside of the boundary layer, a relatively simple general approxi- 
mate solution to the boundary- layer equations has been derived for a 

Prandtl number of unity and for a linear relation between the viscosity 
coefficient and the temperature. This solution not only is convenient 
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to apply but also leads to results which will be sufficiently accurate 
for most practical purposes, including determination of critical Reynolds 
number. 


2. For the particular purpose of determining the location of the 
separation, point in a region of adverse pressure gradient, a relatively 
quick a.nd simple method has been given here, based on the use of seventh- 
degree velocity profiles. The numerical example indicated that the 
results thus obtained should in general be quite accurate. 

3. In the special case of flow near a forward stagnation point, the 
Karman-Pohlhausen method with fourth- degree_ velocity profiles gives 
results of very satisfactory accuracy, even for the critical Reynolds 
number . 

4. The numerical examples given here illustrate the destabilizing 
influences of an adverse pressure gradient and of a Mach number greater 
than zero. 


Polytechnic Institute 
Brooklyn, N. Y., 


of Brooklyn 
May 28, 1951 
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APPENDIX A 

PHYSICAL BOUNDARY -LAYER THICKNESS 


Results in the xt-plane can he transformed into the physical 
xy-plane by means of relation (10). Thus 


y/St 



(Al) 


By inserting velocity profile (l6) and equations (l9a) and (l9b) into 
expression (6a) for (t/T-A, equation (Al) yields: 


X 

St 


= T - 






i + wjh 


u x 2 N 

u m 2 , 


~ t + ft 


f2 Y 3 4 
5/ T 


a 2 


(>-D 


4 a 2 2 5 

T + -5— "T 

5 


- f( 5 + 2a 2 )( 1 - t) t6 + 


f(3 + «a)(l - f) - f =2(5 - 2b 2 ) 
j( 3 + »a) - |(l - ^)(a + | b 2 )]t 8 


(5 + 2a 2 ) 2 - 2a 2 (a + | a 2 ) t 9 - 


|(3 + a 2 )(5 ♦ 2a 2 ) T 10 + 


3j(5 + 2a 2 )(2 + | a a ) + jr(3 + aa)^ 11 - 
| (3 + «a)(2 + | -a)^ + ^( 2 + | =2) 2 ^ 13 - 


(A2) 
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After 5^.(| ) has been obtained, the point g,y in the physical plane 

corresponding to any given point in the mathematical plane can be 

directly determined by means of equation (A2) . 

The boundary-layer thickness S in the physical plane is the value 
of y at which t = 1. Hence, by putting t = 1 in equation (A2), the 
following explicit expression is obtained for the boundary-layer 
thickness: 



I2L -.. A) . 


1 + 


VW 


1 + 


Me/ 1 - 




(o.3950 + 0.021l6a 2 - 0.000622a 2 2 ) 


(A3) 
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APPENDIX B 

DETERMINATION OF SEPARATION POINT 


A seventh-degree -polynomial velocity profile satisfying condi- 
tion (24) in addition to conditions (15) is 


u 

U 1 




lOa^T^ 


(7 + 3a 2 )-7-6 



(Bl) 


where 82 is given by equation (20). With the profile (Bl), the F's 

as defined by equations (12) become (a subscript s is used for the 
present results to distinguish them from the corresponding results based 
on sixth-degree profiles): 


F ls = 0.1156 + 0.00253a 2 - 0.00l454a 2 2 

F 2s = -0.02976a 2 - 0.3125 - Mi 2 (o.428l + 

0.03229a 2 - 0.00l45a 2 2 ) 


(B2) 


F 3s = 1.7500 - 0.5000a 2 


J 


The separation point will be located where 



it 


therefore follows from equation (Bl) that laminar -flow separation will 
occur at the point where = 3.5. By using expressions (B2) and 


assuming, as in the section entitled "General Approximate Solutions," 
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that in the expressions for F]_ s and F 2s ''the quantities a 2 and a 2 ^ 

may, as an approximation, he replaced hy constant values a 2 and a £>2, 
the following approximate solution X s for X(|) of differential equa- 
tion (l4) (with F-l, F 2 , and F3 replaced. by F^, F 2s ,‘ and F 3s , 

respectively) is obtained: 


J 

r s 

0 


3 + _2_ (0 . 0625+0 . 02976a 2 )/_ ■ 

1 ^ t 

7 1 

" F l* 

(0 . 1781+0 . 03229a 2 -0 . 00 l 45 a^ ) 

1 

F ls G | 

ni \b+ j=-{ 0 .O 625 +O.O 2976 a 2 )/^V 

~ ) 18 It ) 

00/ \ 00/ 

r - 1 F is 

[0 . 1781+0 . 03229S2-0. 00145 ^ ) 

1 


Since the main oojecx ox equation (B3) is to determine the separa- 
tion point, the values of a 2 and a 2 ^ may for this purpose be chosen 

as those which a2 and ag^ would actually have at that point. Hence, 
for determining the separation point, the values a 2 = 3*5 and 
a 2 2 = 12.25 are inserted into equation (B3). With 7 = 1.4, expres- 
sion- (25) for Xg(|) then results. 

Since separation occurs where a 2 = 3.5 and since a 2 , which is 

proportional to X, is given by equation (20), the value of X at which 
separation occurs will be that given by equation (26) in the main text. 

Equations (B3) and (25) for X s (l) are valid only under the bound- 
ary condition X 3 = 0 at | = 0. In case the region of adverse pressure 
gradient starts at some point I = ! a other than the leading edge, this 
boundary condition must be replaced by the condition Xg = X a at (• = £ a . 

Solution (25) for Xg(£) must then be replaced by the following solution 
^with a 2 = 3.5 and a 2 ^ = 12.25^: 


32.83 

G 




Xg - 


/u-l \6 . 128 /T-AO . 9366 


. V u 

la ' 00 ' 


IT 


d| + X 


? Ul \7. 128/0^X2. 437 


, u 




^u 1 \7.128/t 1 \2.437 


(B4) 


u 
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The value of is obtained directly from equation ( 22 ) of the main 

text . Thus 


k L 

n 

\3 + ^-(0.08571+0.0190582)/, 

~J V 

\- 24 - --R 0 . 1950 + 0 . 021 l 6 ao- 0 . 000622 a 2 ) 
^ip F l v 2 

T coJ 

F,G 

l- 

\h+ -£-(0.08571+0.019052 2 )/ Tl 

J Fl 

S=l a V 

\ -^(0.1950+0. 02 U 6 a 2 - 0 . 000622 a 2 2 ) 

^= 4 a 


(B5) 


t 



30 


NACA TN 2784 


APPENDIX C 

DETERMINATION OF CRITICAL REYNOLDS NUMBER 


In accordance with the procedure outlined by Lees in reference 11, 
the critical Reynolds number at a given station x for the stability 
of the laminar boundary layer can be determined in the following manner. 
It is required first to determine the value of y/2 such that 


(C1) 

where 



the prime here denotes differentiation with respect to y /l and Z is a 
convenient length.. In the present analysis the quantity 5^ was chosen 

for Z, except in the case of Howarth's series solution (reference 4) of 

the partial differential equations, where the quantity PCyx/ua,)"^^ was 
chosen for Z. 

The appropriate value of y/Z satisfying equation. (Cl) must usually 
be found graphically by trial. After this value has been determined, the 
critical Reynolds number based on l and on local values can be obtained 
from the expression: 
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where the subscript c denotes values at the point where equation (Cl) 
is satisfied. 

For compressible flow it is more convenient to perform calculations 
in the xt-plane than in the xy-plane. With 2 = 8 t , the foregoing rela- 
tions can be expressed directly in the xt-plane by use of the relations: 


5 = dr / a\ 

d(y/l) S(y/Z)\ST/ 


T St 


a 2 

3(y/J) 2 



where 
to t , 


t = t/B^. Using a dot (’) to denote differentiation with respect 
equations (C2) and (C3) become 


a = -jr 
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The quantity y/S^ in terms of t is given hy equation (A2) in 

appendix A) the trial value of. t (instead of y/z) must now he found 
which satisfies equation (Cl). 
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TABLE I 

SEPARATION POINTS CALCULATED BY APIKOXIMAIE METHOD (EQUATIONS (25) AND (26)) 
COMPARED WITH STEWARTSON'S SOLUTION AND WITH OTHER CALCULATIONS 



Method 

5 at separation point for Moo of - 

0 

l 

3 

Method of equations (25) and ( 26 ) 

0.122 

0.113 

0.0768 

Stewartson' s solution, (reference 5) 

.120 

.110 

.0-77 

Method based on sixth-degree profile 




( equation (22)) 




Method based on fourth-degree profile 

.156 


1 



TABLE II 

VALUES OF CRITICAL REYNOLDS NUMBER FOR STABILITY AS CALCULATED 
FROM APPROXIMATE SOLUTION (EQUATION (22)) OF EQUATION (l4) 
COMPARED WITH OTHER SOLUTIONS. FOR ZERO MACH NUMBER 



Method 

(Bj c for of - 

0 

1 

3 

Approximate method, 
sixth-degree profile 
(equation (22)) 

Exact solution of 
equation (l^) 

Hovarth's series 

solution (reference h) 

Flat-plate flow 
(reference 8) 

330 X 10 3 

380 

290 

1500 


isoo(i ♦ |-) 3 

( 1.2 * 01.OJS36 * 0 

( 2.6 ♦ l)(l.l 7i 2 * if 
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TABLE III 

COMPARISON OF CALCULATED SKIN-FRICTION COEFFICIENTS 

[sc - - °l 


Method 

Roo^f /l 

Fourth-degree profile (equation (35)) 

2.390 

Sixth-degree profile (equation (32)) 

2.530 

Solution of Falkner and Skan (reference 14) 

2.468 



TABLE IV 

COMPARISON OF CALCULATED CRITICAL REYNOLDS NUMBERS FOR STABILITY 



Method 

(K„) C S 2 

Fourth-degree profile (equation (35)) 
Sixth-degree profile (equation (32)) 
Solution of Falkner and Skan (reference l4) 

210 X 10 6 
170 
240 
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Figure 2.- Skin- friction coefficients obtained by approximate solution 

fpn iw h^ nn ( OO ^ ^ nf Ard 1 no >nr H ■! f f o>*on+ Hoi q<tiiq + -< nn ( Q/nm + i^vi ^ *1 )i ^ 

\ 1 \ ’ — *- / / V* WJ. U.-UJJMJ. j k-u.uj.uj. ^\-l^UO)UJ.WU ^ J_T f j 

compared -with Howarth 1 b solution for = 0 (reference 4). 

ut/uoo =1-1. 
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Figure 4.- Comparison of fourth- and sixth-degree Telocity profiles 

with solution of Falkner and Skan (reference 14). U]_/u M = 1^, = 0. 
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Figure 6. - Skin- friction coefficients based on approximate solution 
(equation (36)) of ordinary differential equation (equation (l4)); 
fourth-degree velocity profile, up/u^, = £. 
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0 .4 .8 1.2 1.6 2.0 


± r~ L 

L V CTJ W 

Figure 7*- Velocity profiles based on approximate solution (equation ( 36 )) 
of ordinary differential equation (equation (l4)); fourth-degree 
velocity profile, u^ = i; £ = 0.4. 
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